An alternate model for magnetization plateaus in the molecular magnet V15 
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Starting from an antiferromagnetic Heisenberg Hamiltonian for the fifteen spin-f/2 ions in Vis, 
we construct an effective spin Hamiltonian involving eight low-lying states (spin-f/2 and spin-3/2) 
coupled to a phonon bath. We numerically solve the time-dependent Schrodinger equation of this 
system, and obtain the magnetization as a function of temperature in a time-dependent magnetic 
field. The magnetization exhibits unusual patterns of hysteresis and plateaus as the field sweep 
rate and temperature are varied. The observed plateaus are not due to quantum tunneling but 
are a result of thermal averaging. Our results are in good agreement with recent experimental 
observations. 
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PACS numbers: 75.45.+j, 75.60.Ej, 75.50.Xx 

The synthesis of high-nuclearity transition metal clus- 
ters such as Mni2, Feg and V15 jl| has provided an im- 
petus to the study of magnetism on the nanoscale. These 
transition metal clusters are basically isolated transi- 
tion metal complexes involving multi-dentate ligands; the 
chemical pathway between the metal ions in the transi- 
tion metal complex dictates the nature of exchange in- 
teractions. The complex interplay of the topology of ex- 
change interactions, magnetic dipolar interactions and 
spin-lattice coupling has yielded a rich physics on the 
nanoscale which includes quantum tunneling, quantum 
phase interference and quantum coherence [§J3|. Quan- 
tum resonance tunneling is characterized by the observa- 
tion of discrete steps or plateaus in the magnetic hystere- 
sis loops at low temperatures. The signature of quantum 
interference is seen in the variation of the tunnel splitting 
as a function of the azimuthal angle of the transverse field 
for tunneling between M s — — 10 and 10 — n states in 
molecular magnets with ground state spin-10 Q. Quan- 
tum coherence-decoherence studies are important from 
the stand point of application of these systems in quan- 
tum computations Q. 

There have been several models proposed to under- 
stand these phenomena Q. Quantum hysteresis and in- 
terference have largely been studied by using an effective 
spin Hamiltonian with dipolar interactions with a time 
varying external magnetic field. The time evolution of 
the states of the system are carried out within a master- 
equation approach 0. The decoherence phenomena has 
been studied by using a simple two state model in a trans- 
verse magnetic field ^ . Even though most of these clus- 
ters contain a fairly small number of metal ions, the spin 
on the metal ion, at least in the case of Fes and Mni2, 
is fairly large; a full quantum mechanical study of these 
systems is difficult because of the large Fock space dimen- 
sionalities of 1.69 million and 100 million respectively. 
However, the V15 cluster is far more amenable to a rig- 
orous quantum mechanical analysis because of the much 
smaller Fock space (« 33, 000 dimensional) spanned by 



the unpaired spins of the system. A quantitative study of 
these systems requires at least the low-lying states of the 
full spin-Hamiltonian to be evolved in time quantum me- 
chanically as the external magnetic field is ramped with 
time (as is done in experiments). In this letter, we study 
the magnetization of V15 by following its evolution as a 
function of a time-dependent magnetic field at different 
temperatures. The low-lying states are obtained by solv- 
ing the exchange Hamiltonian corresponding to all the 
spins of the system. A spin-phonon interaction is then 
introduced in the Hamiltonian. We thermally average 
the magnetization over the low-lying states after each of 
these states is independently evolved. We find that this 
model reproduces quantitatively all the experimental fea- 
tures found in the magnetization studies of V15 |J. 

The schematic structure of V15 is shown in Fig. 1. 
Structural and related studies on the cluster indicate 
that within each hexagon, there are three alternating ex- 
changes J 800 K which are the strongest in the system, 
and they define the energy scale of the problem. Besides, 
there are weaker exchange interactions between the spins 
involved in the strong exchange and also with the trian- 
gle spins which lie between the hexagons. All the ex- 
change interactions are antiferromagnetic in nature. The 
exchange pathways and their strengths |9| are also shown 
in Fig. 1. What is significant in the cluster is the fact 
that the spins in the triangle do not experience direct 
exchange interactions of any significance. The exchange 
Hamiltonian of the cluster is solved using a valence bond 
basis in each of the total spin subspaces, for all the eigen- 
values. It is found that two spin-1/2 states and a spin-3/2 
state are split-off from the rest of the spectrum by a gap 
of 0.6 J 0. These eight states almost exclusively corre- 
spond to the triangle spins and they are the only states 
which will make significant contributions to sub-kelvin 
properties. We therefore set up an effective Hamiltonian 
in the Fock space of the three spins. We find that the 
form 
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+ 6 Sx ■ S 2 x S 3 , (1) 

where e = -4.781808, a = -0.001491, and S = 0.015144 
in units of the exchange J (see Fig. 1), reproduces the 
low-lying eigenstates to numerical accuracy. Note the 
unusually large value of the three-spin interaction; this 
term has not been considered in the earlier literature on 
this subject (such as Ref. [9]), and it is essential for the 
spin-3/2 state to lie between the two spin-1/2 states as 
is found to be the case in this system. 

The direct spin-spin interaction terms permitted by 
the C3 symmetry are given by 

H dip = 7 [ (sl + S 3 _) + i (S 3 + -S 3 _)}. (2) 

We have also introduced a coupling between the spin 
states of the cluster and the phonons. The spin-phonon 
interaction Hamiltonian which preserves the C3 symme- 
try is phenomenologically given by |Tcj 

H sp - ph = q(b + tf)[ (S 2 + + S 2 _) + i(Sl - S 2 _) 

+ (S 2 Z ~ \S 2 ) } , (3) 

where q is the spin-phonon coupling constant, b (b*) is 
the phonon annihilation (creation) operator, and Huj is 
the phonon frequency. For simplicity, we have assumed 
a single phonon mode although the molecule has various 
possible vibrational modes. The form of the interaction 
in Eq. (p|) means that the phonons couple only to states 
with spin-3/2. We have restricted the dimensionality of 
the Fock space of the phonons to 15 considering the low 
temperatures of interest. 

The evolution of the magnetization as a function of 
the magnetic field has been studied by using the total 
Hamiltonian H to tai, given by 

Htotai — H sp - sp + Hdip + Hsp-ph + hcj(b'b + 1/2) 
+ h z (t)S z + h x (t)S x , (4) 

where we have assumed that besides an axial field h z (t), 
a small transverse field h x (t) could also be present to 
account for any mismatch between the crystalline z-axis 
and the molecular z-axis. The numerical method involves 
setting up the Hamiltonian matrix in the product basis 
of the spin and phonon states \i,j >, where \i > cor- 
responds to one of the eight spin configurations of the 
three spins, and j varies from to 14, corresponding to 
the fifteen phonon states retained in the problem. The 
values we have assigned to the different parameters are 
7 = 10~ 3 , q = 10~ 4 and HbJ = 1.25 x 10~ 4 , all in units of 
the exchange J (see Fig. 1). 

To study the magnetization phenomena, we start with 
the direct product eigenstates of H sp - sp and Tiu>(Wb + 
1/2), and independently evolve each of the 120 states ipij 
by using the time evolution operator 

tp(t + At) = e- iHt ° talM ' h il)(t) . (5) 



The evolution is carried out in small time steps by ap- 
plying the evolution operator to the state arrived at in 
the previous step. The magnetic field is changed step- 
wise in units of 0.015T. At each value of the magnetic 
field, the system is allowed to evolve for 300 time steps 
of size At, before the field is changed to the next value. 
At every time step, the average magnetization < M(t) > 
is calculated as 

JL J± e -/3K+M*)TO<] e -Phu>v+i/2) 

< M (t)> = y: E — — z- h — 

1 ^ spin ^ph 

<^M*)l&liM*) > , (6) 

where Wi and mi are the eigenvalues and magnetizations 
of eigenstates of H sp _ sp , f3 — l/fc^T, Z p h is the phonon 
partition function of Tiuj(b^b+l/2), and Z sp i n is the parti- 
tion function of H sp ^ sp in the presence of the axial mag- 
netic field. In Fig. 2, we show the energy level ordering of 
the effective spin Hamiltonian and the effect of the mag- 
netic field on the eigenvalue spectrum. We also show the 
couplings between various states brought about by the 
magnetic dipolar terms and the spin-phonon terms; note 
that the spin-1/2 and spin-3/2 states are not connected 
to each other by these terms. 

In Fig. 3, we show the hysteresis plots of the system for 
different temperatures. We see that at low temperatures, 
the plateaus in the hysteresis plots are very pronounced. 
The plateau width at < S z >= —0.5 corresponds to 2.64 
T which is in excellent agreement with the experimental 
value of 2.8 T § assuming that J = 800 K. We also 
find that the plateau vanishes above a temperature of 
0.8 K which is also in excellent agreement with the ex- 
perimental value of 0.9 K [Q. The inset in figure shows 
the temperature variation of the plateau width. We note 
that the plateau width falls off rapidly with temperature, 
and an exponential fit to W — A exp(— T/fi) (see Fig. 
3) gives the characteristic temperature O to be 0.2 K. 
This small value of f2 is because there are no large bar- 
riers between the different magnetization states in this 
system, unlike the high spin molecular magnets such as 
Mn 12 §. 

We also observe that when the field is swept more 
rapidly, there are additional plateaus at intermediate val- 
ues of magnetization. For example, in Fig. 4 the field 
sweep rate is increased by a factor of five compared to 
Fig. 3, and we find a small plateau of width 0.03 T 
near H = 0.15 T at a value of < S z >= -0.375. This 
is because near that field, some of the spin-3/2 states 
become degenerate in energy; subsequently, as the mag- 
netic field is increased, the system stays locked in some 
of those states if the sweep rate is too high. This plateau 
vanishes on warming the system slightly. We have also 
studied the effect of cycling the field on the plateaus. In 
Fig. 5, we show the magnetization as a function of field 
at different temperatures and sweep rates. It is interest- 
ing to see that the system does not have either remnance 
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or coercivity; all the hysteresis plots pass through the ori- 
gin. The effect of varying the rate of scanning the field is 
also shown in Fig. 5. We find that as the scanning rate 
increases, the hysteresis in the plot of magnetization vs 
field decreases, and the plateau feature is almost identical 
in both the scanning directions. This could be because of 
the slow relaxation of the magnetization which is indeed 
the reason why the plateaus occur in the first place. We 
also find that the transverse field term does not affect 
any of our results significantly. 

To summarize, we have derived an effective Hamilto- 
nian from the exchange Hamiltonian of the full V15 sys- 
tem. In the presence of a time varying magnetic field, the 
states of the effective Hamiltonian are allowed to evolve 
under the influence of magnetic dipolar interactions and 
a spin-phonon coupling. During the time evolution, the 
magnetization is followed as a function of the applied 
magnetic field. The calculated M vs H plots show mag- 
netization plateaus at low temperatures. The width of 
the plateau at low temperature as well as the tempera- 
ture at which the plateau vanishes are in excellent agree- 
ment with experimental values. It is also shown that the 
number of plateaus observed depends upon the scanning 
speed of the magnetic field. When the magnetic field is 
cycled, the hysteresis plots pass through the origin indi- 
cating the absence of remnance and coercion. The hys- 
teresis is pronounced for slow scanning speeds. From our 
results, it appears that the magnetization plateaus in V15 
is not a consequence of quantum resonant tunneling but 
is a result of thermal averaging. We also find that the 
magnetization does not show any oscillation with time 
during evolution indicating the absence of quantum tun- 
neling. 
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J- 800K; J ~ 150 K; 1 ~ 300 K 

FIG. 1. Schematic exchange interactions in a V15 cluster. 
There is no direct exchange interaction amongst the triangle 
spins. Interactions not shown explicitly can be generated from 
the C3 symmetry of the system. 
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FIG. 2. (a) Eigenstates of the effective spin Hamiltonian 
H S p- S p, (b) Eigenstates in the presence of a moderate ax- 
ial field. Arrows show the states connected by the dipolar 
terms and the transverse field, (c) is the same as (b) but in 
a stronger field, (d) describes the effect of spin-phonon terms 
(shown by arrows with broken lines) on (c). 



FIG. 4. Magnetization as a function of axial field for a 
faster sweep rate at three different temperatures. 




FIG. 3. Plot of magnetization vs axial field at different tem- 
peratures. Inset shows plateau width as a function of temper- 
ature (full circles). Triangles in the inset correspond to the 
values from the fit to an exponential function. 



FIG. 5. Magnetization vs axial field for a full cycling of the 
field at different temperatures and sweep rates. 
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